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The thermodynamic quantities (spin-spin correlation lunctions (SoS n ), correlation length f, 
spin susceptibility Xj an d specific heat Cv) of the frustrated one-dimensional Jx-Ji Heisenberg 
ferromagnet with arbitrary spin quantum number S below the quantum critical point, i.e. for 
Ji < I -A 1/4, are calculated using a rotation-invariant Green- function formalism and full diago- 
nalization as well as a finite-temperature Lanczos technique for finite chains of up to N — 18 
sites. The low-temperature behavior of the susceptibility x an d the correlation length £ is well 
described by X = f^ 4 (| Ji| - 4J 2 ) T~ 2 + AS 5/2 (| Ji| - 4J 2 ) 1/2 T~ 3/2 and £ = S 2 (|Ji| - 4J 2 ) T' 1 + 

BS 1/2 (\Ji\- 4J 2 ) 1/2 T~ 1/2 with A « 1.1 . . . 1.2 and B « 0.84 . . . 0.89. The vanishing of the factors 
in front of the temperature at J 2 = |Ji|/4 indicates a change of the critical behavior of \ and £ at 
T — ► 0. The specific heat may exhibit an additional frustration-induced low-temperature maximum 
when approaching the quantum critical point. This maximum appears for S =1/2 and S = 1, but 
was not found for S 1 > 1. 



PACS numbers: 

I. INTRODUCTION 

The one-dimensional (ID) Heisenberg model with fer- 
romagnetic nearest-neighbor (NN) exchange J\ < 
and frustrating antiferromagnetic next nearest-neighbor 
(NN) exchange </ 2 > has recently attracted much 
attention, see, e.g., Rcfs. 

[Ml 

This model may 

serve as the minimal model to describe the phys- 
ical properties of edge-shared chain cuprates, e.g. 
Ca 2 Y 2 Cu 5 Oio, Li 2 ZrCu0 4 , and Li 2 Cu0 2 ^^ On the 
other hand, several materials considered as ID ferro- 
magnets, such as TMCuC[(CH 3 ) 4 NCuCl 3 ] A CsNiF 3 ^2 
Ni(4-cyanopyridine) 2 Cl2^£ or deposited Co chains^ 
might have also a weak frustrating NNN exchange J 2 < 
— Ji/4. Moreover, recent investigations suggest that 
Li 2 Cu0 2 as well as Ca 2 Y 2 Cu50io are quasi-lD spin sys- 
tems with a dominant ferromagnetic J\ and weak frus- 
trating antiferromagnetic J 2 so that the inchain spin-spin 
correlations are predominantly ferromagnetic! 12 ' 13 ' 15 

The corresponding Hamiltonian of the ID Ji-J 2 
Heisenberg model considered in this paper is given by 



H 



(JlSjSj+l + J 2 SiS; + 2 



(1) 



where i runs over all lattice sites and S 2 = S(S + 1). 
We set Ji = — 1 and consider J 2 > 0. Although some of 
the above mentioned materials, namely CsNiF3fii Ni(4- 
cyanopyridinc^Cl^ or deposited Co chains^, represent 



ID ferromagncts with spin quantum number S > 1/2, so 
far the focus of recent theoretical studies has been on 
the S — 1/2 case. Since Haldanc's famous paper— we 
know that in ID Heisenberg systems the spin quantum 
number may play a crucial role. Recently it has been 
found that the frustrated model |T]) in the extreme quan- 
tum case (5=1/2 may exhibit a different behavior near 
the quantum critical point than the model for S > 1/2, 4 - 
Moreover, in Ref.[2l|it has been found that for the unfrus- 
trated ID quantum ferromagncts a characteristic field- 
induced low-temperature maximum in the specific heat 
exists only for the small spin qu antum numbers (5=1/2 
and S = 1, see also Refs. QjlalU 

In the present paper we discuss the thermodynam- 
ics of the model ([1} for arbitrary spin quantum num- 
ber S and focus on the parameter regime J 2 < |Ji|/4, 
where the ferromagnetic ground state is realized. It 
has been recently demonstrated for low-dimensional frus- 
trated ferromagncts 1 7 i 1 1 i 24 i 25 that the frustrating J 2 may 
influence the thermodynamics substantially. In the ID 
system a change in the low-temperature behavior of the 
susceptibility and the correlation length as well as an ad- 
ditional low-temperature maximum in the specific heat 
have been found when approaching the zero-temperature 
critical point at J 2 = | J\ |/4 from the ferromagnetic side»£ 
In particular, a different critical behavior of the suscep- 
tibility and the correlation length has been found for 
J 2 < |^i |/4 and at J 2 = \J x \j\^ 
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As in our previous papers on frustrated spin- 
1/2 ferromagnets ; 7 i 1 1 1 24 we use a second-order Green- 
function technique 2 ^ to study the influence of the spin 
quantum number S on the thermodynamic properties 
of the model ([l). This technique has been applied 
successfully to several low-dimensional quantum spin 
systemsiLlUlraiJHtt 

To extend our Green's function theory! to S > 1/2, we 
follow Ref. where the Green's function technique was 
applied to a ferro- and antiferromagnetic layered square 
lattice with arbitrary S. We complement the Green's 
function results by full exact diagonalization (ED) and 
finite-temperature Lanczos (FTL) technique data for fi- 
nite systems of up to N = 18 lattice sites. 



II. FULL DIAGONALIZATION AND 
FINITE-TEMPERATURE LANCZOS 
TECHNIQUE FOR FINITE LATTICES 

Using Schulcnburg's spinpaclt^ and exploiting the lat- 
tice symmetries and the fact that S z = J2i commutes 
with H we are able to calculate the exact thermodynam- 
ics for periodic chains of up to N = 14 spins for spin 
quantum number S = 1. For N = 8 sites we found the 
exact thermodynamics up to S = 2. Clearly, the full 
ED of finite systems for S > 1/2 is less efficient as for 
S = 1/2, since the accessible system size N decreases 
with increasing of S. 

In addition to the full ED we have also applied a FTL 
technique, see e.g. Refs. I33ll34l . This method allows an 
accurate calculation of thermodynamic quantities down 
to low temperatures for S = 1 up to N = 18. 



(S7S+)S+ + \(S+S7)S+ for products with two coincid- 
ing sites which occur for S > L 21 ' 27 i 29 ' 30 Note that for 
J 2 < — Ji/4, where the ground-state is ferromagnetic, 
the vertex parameters a and A can be assumed in a good 
approximation to be independent of the range of the as- 
sociated spin correlators, cf. Ref. 0. Then we obtain 



-5+ = lo 2 S+ and 



Q 1 — a 1 1 tu 9 9' 



(3) 



with 



lo 2 = ^2 JnJ m (1 - cos nq) [K n<m + AaC n (1 - cosmg)] , 

n.m(=l,2) 

(4) 

where K n . n = |5(5 + 1) + 2\C n + 2a(C 2n -3C n ), 
K 1)2 = 2a (C 3 ~ d), and K 2A = K\ 2 + 4a (d - C 2 ). 
From the Green function Q the correlation functions 
= W C q e lqn are determined by the spectral 
theorem^ 2 - 



M„ 



C g = (S+SZ q ) = ^-[l + 2n(uj q )] 



(5) 



where n{u) q ) = (e Wq ' — l) is the Bose function. Using 
the operator identity Sf = SfSr - Sf + (Sf ) 2 , we get 
the sum rule C = ^ C i = jS(S+l). Following 
Ref. l30l . as an additional equation to determine the vertex 
parameters we consider the ratio 



r(T) 



A(T) - A(oo) 
a(T) — a(oo) 



r(0) 



(6) 



III. SPIN-ROTATION-INVARIANT 
GREEN-FUNCTION THEORY 

To study the thermodynamics of the model ([1]) for 
arbitrary N we use the spin-rotation-invariant Green- 
function method (RGM) ZZ^ZkMgo The relevant ther- 
modynamic quantities can be determined by calculating 
the two-time commutator Green function^ ((<S+; SZ q ))u, 
which is related to the transverse spin susceptibility by 
X q ~ (uj) = —((S q ;SZ q ))u ■ Using the equations of mo- 
tion up to the second step and supposing rotational 
symmetry with (Sf) = 0, we obtain uj 2 ((S+; SZ g ))u = 
M q + ((-S+-S 
-S, 



'q L I q 

expression 



_,))„ with M q = ([[S+,H],SZ q ]) and 
For the moment M q the exact 



[[S+,H\ ,H\ 



Mr, 



E 



J n C n (1 - cosnq) 



(2) 



holds, where C„ = {S£S~) = 2(S$S*}. The sec- 
ond derivative — S q is approximated in the spirit of 
Refs. HH-ilil-lil i.e., in -S+ we use the decoupling 
S+S+S- = a(S+S-)St + a(S+S~)S+, and S+S7S+ = 



as temperature independent, where A(oo) = 1 — 4 g(f +1 ) 
and a(oo) = 1, see Ref. [29|. The uniform static spin 
susceptibility \ — li m <?^o Xqi where \q — Xq ( w — 0) an d 
X q (w) = \xt~ M, is given by 



X 



n 2 J n C n ; A 



71=1,2 



n,m(=l,2) 



(7) 

The correlation length £ for a ferromagnet can be cal- 
culated from the expansion of the static spin suscepti- 
bility around q = (see, e.g., Refs. [26| and liHI). Xq = 
xj (l + £, 2 q 2 )- The ferromagnetic long-range order, oc- 
curring in the ID model at zero temperature only, is 
related to the condensation term C (see Ref. [26h via 
C n (0) = £ £ ?(/0) (M q /2uj q ) e^ + C. Equating this ex- 
pression to the exact result C n (0) = ^S6 ni o + ^S 2 leads 
to C(0) = |S* 2 and M 9 (0)/2w g (0) = This requires 

a(0) = f and K n , m {0) = [cf. Eqs. © and Q] which 
leads to A(0) = 2 — ^. The parameter A in Eq. ([7]) is 
zero at T = 0, i.e., x diverges as T according to the 
ferromagnetic phase transition. Moreover, using these 
results we find at zero temperature uj q = 2p s q 2 (\q\ <C 1), 



Ps 



S 



(| Ji| — 4J2), where p s is the spin stiffness. 
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Since we will compare the RGM data with ED data for 
finite lattices, see Sec. IIV1 we have to adopt the RGM to 
finite A. In this case the quantity C is not related to the 
magnetisation and stays nonzero in the whole tempera- 
ture region. In Ref. [29| it was shown that C = 2T\/N. 

Finally, to solve the equation of motion and to evaluate 
the thermodynamic quantities we have to determine the 
correlators C; (i = 1, . . . ,4) and the vertex parameters a 
and A (for finite systems, also C) as the numerical solu- 
tion of a coupled system of six (seven) non-linear alge- 
braic self-consistency equations: C; according to Eq. §5§ 
including the sum rule Co = ^, and the ratio Eq. (fr?]). 
To find the numerical solution of the RGM equations 
for T > 0, we start at high temperatures and decrease 
T in small steps. Below a certain (low) temperature 
Tq(J2,S) no solutions of the RGM equations (except at 
T = 0) could be found, since the quantity A(T, J2, S) in 
Eq. (J7J becomes exponentially small which leads to nu- 
merical uncertainties. We find that To (J 2 , S)/S(S + 1) 
is of the order of 10~ 3 for all values of J 2 and S con- 
sidered here. Moreover, we find that the RGM solution 
at very low temperatures becomes less trustworthy for 
J 2 approaching the quantum critical point J 2 = |Ji|/4, 
cf. also Refs. andlTJ. Therefore, below we will present 
RGM results for J 2 < 0.22| JJ only. 



IV. RESULTS 

Motivated by the high-temperature behavior of the 
physical quantitie o 37 ! 38 we use as the renormalized tem- 
perature scale t = T/S(S + 1), since for large temper- 
atures physical quantities exhibit ^-dependence which is 
independent of S. 



A. Spin-spin correlation functions 

Let us first consider the spin-spin correlation functions 
(S S„) depicted in Fig. Q] for NN, NNN and tenth- nearest 
neighbors. Analogously to the frustrated 5 = 1/2 chain^ 
with increasing frustration the correlation functions de- 
crease more rapidly. Obviously the spin quantum number 
S has only a small influence on the correlation functions 
as functions of t. Interestingly, an increase of S yields a 
slight weakening of the short-range spin-spin correlation 
at fixed t > (see upper and middle panels in Fig. [J), but 
an enhancement of larger-distant correlations (see lower 
panel in Fig. [TJ. Moreover, the NNN correlation function 
for J 2 = 0.2 becomes negative at t = t ps 1.35, 1.37, 1.38 
for S = 1,3/2,2, respectively. (Note that for 5=1/2 
the corresponding temperature is to ~ 1.27^) 



B. Susceptibility and correlation length 

The behavior of the correlation functions is reflected in 
the susceptibility \ and the correlation length £ shown in 
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FIG. 1: Spin-spin correlation functions (S Si)/S 2 (NN), 
(S S 2 )/S 2 (NNN), and (S Si )/5 2 obtained by the RGM 
for spin quantum numbers S — 1 (solid), 3/2 (long-dashed), 
and S — 2 (short-dashed) and frustration parameters J2 = 
(black), 0.1 (red), and 0.2 (blue). 



Figs. [2] and [3l We illustrate the influence of frustration 
for a fixed spin quantum number S = 1, whereas in the 
insets we show the influence of S for a fixed frustration 
parameter J 2 = 0.1. Since the ground state is ferromag- 
netic, both quantities diverge at T = 0. With increasing 
of frustration the rapid increase in x and £; is shifted to 
lower temperatures. The comparison of RGM and ED 
results for finite chains of A = 12 presented in Fig. [5] 
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FIG. 2: Susceptibility \ f° r spin quantum number S = 1 and 
frustration parameters J2 = (black), 0.1 (red), 0.2 (blue), 
and 0.22 (violet) [solid lines: RGM for TV — > 00, dashed lines: 
RGM for N = 12, symbols: ED for N = 12]. Inset: Suscepti- 
bility x for J 2 = 0.1 and 5 = 1/2, 1, 3/2, and 2 (from bottom 
to top). The data for 5 = 1/2 are taken from Ref. 0. 
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FIG. 3: Correlation length £ for spin quantum number 5 = 
1 and frustration parameters J2 = (black), 0.1 (red), 0.2 
(blue), and 0.22 (violet). Inset: Correlation length £ for J2 = 
0.1 and 5 = 1/2, 1, 3/2, and 2 (from bottom to top). The 
data for 5 = 1/2 are taken from Ref. 0. 



demonstrates a very good agreement of the susceptibility 
data obtained by both methods. It can also be seen that 
the renormalizcd temperature £, where finite-size effects 
become relevant, is shifted to lower values of t with in- 
creasing frustration (in Fig. [3]the curves for N = 12 and 
TV — > 00 for J 2 = 0.2 and J 2 = 0.22 almost coincide). 
The curves shown in the insets again illustrate that the 
influence of the spin quantum number is small; only for 
5=1/2 the curves are noticeably separated from the 
bundle of curves for 5 > 1/2. Moreover, it is evident 
that x an d £ at fixed t and J 2 increase with growing 5. 
Next we consider the critical behavior of x an d £ for 
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FIG. 4: RGM results for fx? 2 / 'S 4 (upper panel) and £T/S 2 
(lower panel) versus VT for spin quantum numbers 5 = 1 
(solid), 3/2 (long-dashed), and 5 = 2 (short-dashed) and frus- 
tration parameters J2 = (black), 0.1 (red), and 0.2 (blue). 



T^Oby analyzing the RGM data at low temperatures. 
From Ref. [13 it is known that within the RGM the low- 
temperature behavior of x and £ of the unfrustratcd fer- 
romagnet (J2 = 0) is given by limT^o X^ 2 = 25 4 |Ji|/3 
and liniT^o £T — S 2 \Ji\. These results agree with those 
obtained by modified spin-wave theor y 39 ' 40 Moreover, 
for 5=1/2 they coincide with the exact Bethe-ansatz 
analysis 4ML (Note that x defined in Refs. [H and EH is 
larger by a factor of four than x given by Eq. 0.) Includ- 
ing higher-order terms in T for ID Hcisenberg fcrromag- 
nets the low-temperature behavior of the susceptibility 
and the correlation length reads \T 2 = yo + TJiVt + 
y 2 T + C(T 3 / 2 ) and fT = x + x x VT + x 2 T + C(T 3 / 2 ), 
see, e.g., Refs. [7l l39l44l[ The dependence of %T 2 and 
£T on T 1 / 2 is shown in Fig. 2] for various values of 5 
and J2. The validity of linear relations xT 2 oc T 1 / 2 and 
£T oc T 1 / 2 at low temperatures is clearly seen. To de- 
termine the parameters yo, y\, y 2 , xq, x\, and x 2 , we 
follow the lines of Ref. and fit the RGM results for x 
and £ using the relations given above. For the fits we use 
low-temperature data points between Tq and To + T cut , 
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where To is the lowest temperature for which solutions 
of the RGM equations can be found, see Sec. IIIIl For 
T C ut we chose T cut = 0.005, cf. Ref. 0. By analyzing 
data for J 2 = and 5 = 1, 3/2, 2 and 5 we have numer- 
ically confirmed the relations yo = 25 4 /3 and xo = S 2 
with an accuracy of at least three digits. For the next-to- 
lcading coefficients we adopt the relations yi/S 4 = aoS ai 
and Xi/S 2 = faS^ 1 found by modified spin- wave the- 
ory for the unfrustrated modeL 39 i 40 We find ao = 1.106, 
ai = -1.497, fa = 0.834, and fa = -1.495, which is 
in good agreement with the results of Refs. [H and Eol . 
where a = 0.824, a x = -1.5, fa = 0.412, and fa = -1.5 
were reported. 

Now we determine the parameters yo-, yi, y 2 , %0y 
and X2 for finite frustration J2 > 0. The results for 
the leading coefficients yo and Xo are shown in the in- 
sets of Fig. [5] Both coefficients y and xo obey with 
high precision the linear relations yo/S 4 = | (| J\ \ — 4J 2 ) 
and Xo/S 2 = \J\ \ — 4J 2 , i.e., the leading coefficients de- 
crease with growing frustration and vanish finally at the 
quantum critical point J| = |Ji|/4, thus indicating a 
change of the critical exponent of x an d £ at J|. This 
linear decrease of yo and Xq found by fitting the low- 
temperature behavior of x an d £ is t ne same as that 
obtained analytically for the zero-temperature spin stiff- 
ness p s , see Sec. IIIIl This correspondence between the 
spin stiffness and the divergence of the susceptibility 
and the correlation length is in accordance with general 
argument o 1 concerning the low-temperature physics 
of low-dimensional Hciscnberg fcrromagncts. Moreover, 
these linear relations agree with the findings for 5 = 1/2 
(Ref. 0) and 5 -> 00 (Ref. M). 

Next we determine the coefficients y\ and as func- 
tions of J2 again using the fit of the RGM data for 5 = 1, 
3/2,2 and 5 described above. Moreover, wc reanalyze our 
data for 5 = 1/2 from Ref. 0. Using an improved com- 
puter code we are able to add 5 = 1/2 data for J 2 = 0.22 
not given in Ref. 0. Similar as for yo/S 4 and xo/S 2 , 
the data for y\ and x\ shown in Fig. [5] yield evidence 
for a universal J2 dependence of yi/5 5 / 2 and xi/S 1 ^ 2 . 
However, the coincidence of the data points is less pro- 
nounced (in particular for 5 = 5) than for yo and Xq- 
The dependence of yi/5 5 / 2 and x\/S x l 2 on J 2 is clearly 
not linear. As suggested by the behavior of yi/5 5 / 2 and 
xijS 1 ! 2 at larger J 2 we chose /(J2) = a yj\ — bJ 2 to fit 
the data points presented in Fig. [S] These fits suggest, 
that yi and x\ also vanish approaching the quantum crit- 
ical point. Indeed, we find for the fit parameter b ~ 4.08 
for most of the fits; only for 5 = 1/2 we have b w 4.20. 
For comparison, we also show the previous quadratic fit 
for 5 = 1/2 of Ref. 0, where the data point at J 2 =0.22 
was not included. It is obvious that the square-root fit is 
more reasonable than the qadratic fit (which leads to a 
finite yi and x\ at J2 = | Ji|/4). Moreover, a square- root 
dependence of y\ and X\ on the exchange couplings is also 
suggested by modified spin- wave theory 39 i 40 i 44 Based on 
the results for y , yi, Xq, and X\ discussed above wc fi- 
nally argue that the low-temperature behavior of the sus- 




FIG. 5: Next-to-leading coefficients y\/S 5 ^ 2 (upper panel) 
and Xx/S 1 / 2 (lower panel) in dependence on J2 obtained by 
fitting the low-temperature data of \ an d £ (see main text). 
The symbols represent the data points for S = 1/2, 1, 3/2, 
2, and S = 5, and the lines of the same color show the cor- 
responding fit curves of these points using the fit function 
/(J2) = ay/1 — &J2. The thick black solid line shows the 
quadratic fit for 5 = 1/2 without the data point at J2 = 0.22 
as used in Ref. 0- 

Insets: Leading coefficients yo (upper panel) and £0 (lower 
panel) in dependence on J2 obtained by fitting the low- 
temperature data of x an d £ (see main text). The values for 
yo/S 4 as well as for xo/S 2 practically coincide for S = 1/2, 
1, 3/2, 2, and 5 = 5 used to determine yo/S 4 and Xo/S 2 . 



ceptibility and the correlation length for J 2 < |»/i|/4 is 
well described by 

XT 2 = ^5 4 (|Ji| - AJ 2 ) + AS 5/2 y / {\J 1 \ ~4J 2 )VT (8) 
and 

£T = 5 2 (\fa\ - 4J 2 ) + BS^ 2 ^(\,h\ - 4J 2 )VT (9) 

with A w 1.1 (A fn 1.2) and B w 0.84 (B w 0.89) for 
5 = 1/2, 1, 3/2, and 2 (5 = 5). Therefore, we conclude 
that, in accordance with previous results for 5=1/2 and 
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S = oo, sec Refs. @ and[25l. the critical behavior of x an d 
£ is changed at the quantum critical point over the entire 
range of spin quantum numbers, 1/2 < S < oo. We 
may speculate, that the expansion \T 2 = yo + y\ 
y 2 T + 0(T 3 / 2 ) and £T = x + x x Vf + x 2 T + C(T 3 / 2 ) 
valid for J 2 < | «/i | / 4 breaks down at the quantum critical 
point J 2 = \Ji |/4, and the critical behavior found for 
the classical system at Jf, x ^ T -4 / 3 and £ oc T -1 / 3 , 
may hold also for the quantum model. However, recent 
numerical studies^ 4 - indicate a slightly different critical 
behavior for 5 = 1/2. 

C. Specific heat 

Finally we discuss the specific heat Cy, see Figs. |H] and 
[7J With increasing of 5 there is a shift of the broad maxi- 




T/S(S+1) 

FIG. 6: Specific heat for 5=1/2 (green), 5 = 1 (black), 
5 = 3/2 (red), and 5 = 2 (blue) obtained by RGM for J 2 = 
(solid), J2 = 0.2 (long-dashed), and J2 = 0.22 (short-dashed). 
Inset: ED results for N = 8 and 5 = 1 (black), 3/2 (red) and 
2 (blue) for J 2 = 0.2 (solid) and ,h = 0.22 (dashed). 



mum typical for spin systems to lower values of the renor- 
malized temperature t = T/S(S + 1) and an increase of 
its height, see Fig. [5] For fixed 5 the frustration J 2 leads 
to a decrease of this maximum. Of particular interest 
is the low-temperature behavior of the specific heat. In 
Ref. an additional frustration- induced low-temperature 
maximum in the specific heat was found for the 5 = 1/2 
model studied by RGM and ED. From Fig. [5] it is obvi- 
ous that for higher values of 5 this additional maximum 
disappears. Note that this observation is similar to the 
findings in Ref. [2l| for a field-induced low-temperature 
maximum in Cy , which appears for low values of 5 only. 
To be more specific, we find that for 5 > 1 no extra maxi- 
mum appears in the RGM data for Cy in the whole range 
of < J 2 < 0.22 accessible within the RGM approach. 
The frustration-induced low-temperature maximum for 
5 = 1/2 appears in the RGM for J 2 > 0.16 (Ref. 0). For 
5=1 we find that a double-maximum structure in Cy 
calculated by RGM appears for J 2 > 0.21. 

These RGM based results are supported by finite-size 
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FIG. 7: Exact diagoalization (ED) and finite-temperature 
Lanczos (FTL) results for the low-temperature maximum in 
the specific heat for 5=1 and J2 = 0.2 (upper panel) and 
,h = 0.22 (lower panel) for chain lengths N = 10, 12, 14, 16 
and 18. For TV = 14 we present both the exact ED and the 
approximate FTL data to demonstrate the accruracy of the 
FTL technique down to very low temperatures. 



data (see inset of Fig. [6]), i.e. we find also in the ED 
data for N = 8 a clear tendency to suppress the extra 
low-temperature maximum if 5 increases. In more detail 
we analyze finite-size data for 5=1, where the larger 
systems are accessible by numerical methods. In Fig. [7J 
we illustrate the finite size-effects in Cy for 5 = 1 and 
J 2 = 0.2 and 0.22. It is obvious that the extra low- 
tcmpcrature maximum is appreciably affected by finite- 
size effects. However, from Fig. [7J it is also evident 
that the height of the extra low-temperature maximum 
as well as the nearby minimum behave monotonously 
with N. Hence a finite-size extrapolation of the height 
C m ax{N) of the extra maximum and the related min- 
imum C m i n (N) is reasonable. Analogously to Ref. 
we find that a(N) = a,Q + a\/N 2 + a 2 /N 4 as a reason- 
able extrapolation scheme. The results of the extrap- 
olation are shown in Fig. [JJ by black [for the extrapo- 
lated value C max {po)\ and red [for the extrapolated value 
of C m in (00)] filled circles. These extrapolated data for 
Cmax and C m i n indicate, that there is most likely no ex- 
tra low-temperature maximum for J 2 = 0.2, but such a 
maximum appears for J 2 = 0.22. Hence, the finite-size 
data support the RGM predictions that the double-peak 
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structure appears for S = 1 at J 2 > 0.21. 

V. SUMMARY 

In this paper we have studied the influence of a frus- 
trating NNN coupling J 2 on the thermodynamics of the 
ID spin-S* J1-J2 Heisenberg model with ferromagnetic 
Ji and antiferromagnctic J 2 with J 2 < |Ji|/4. For that 
we have used the RGM for infinite chains and full ED 
as well as the FTL technique for finite chains. We find 
a universal dependence of the thermodynamic quantities 
on the renormalized temperature t = T/S(S+1) at large 
temperatures. At low temperatures such a universal be- 
havior is found for the critical properties of the suscep- 
tibility xT 2 /S 4 = (2/3)(| Ji| - 4J 2 ) and the correlation 



length £T / S 2 = \Ji\ — 4J 2 , he., the critical exponents 
of x an d £ for T — > are not changed by a frustrating 
J 2 < 0.25|Ji|. However, our data suggest that at the 
quantum critical point J 2 = |</i|/4 the critical behavior 
of x an d £ is changed. 

For S = 1/2 and S — 1 an additional low-temperature 
maximum in the specific heat Cy emerges when J 2 
approaches the quantum critical point. Since we did not 
observe such an additional maximum in Cy for 5 > 1, 
it can be attributed to strong quantum fluctuations 
present at small values of S. 
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